Quantum DifFusion and Localization in Disordered Electronic Systems 
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The diffusion of electronic wave packets in one-dimensional systems with on-site, binary disorder 
is numerically investigated within the framework of a single-band tight-binding model. Fractal 
properties are incorporated by assuming that the distribution of distances £ between consecutive 
impurities obeys a power law, P{£) ~ . For suitable ranges of q, one finds system- wide anomalous 
diffusion. Asymmetric diffusion effects are introduced through the application of an external electric 
field, leading to results similar to those observed in the case of photogenerated electron-hole plasmas 
in tilted InP/InGaAs/InP quantum wells. 

PACS numbers: 72.15.Rn,71.23.-k,71.23.An 
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I. INTRODUCTION 

The work reported in this paper is motivated by 
the experimental observation of asymmetric diffusion 
of electron-hole plasmas in semiconductor quantum 
wellsi*2,»2,, as well as by attempts at a theoretical expla- 
nation of such results, on the basis of models for fractal 
diffusionii^. In Ref. [l|, electron-hole plasmas were photo- 
generated in intrinsic InP/InGaAs/InP single quantum 
wells, and their diffusion was observed through measure- 
ments of photoluminescence intensity profiles. Two types 
of heterostructures were examined: the first was grown in 
a direction normal to the [001] crystallographic axis, and 
the second was grown in a direction tilted by 2° relative to 
that of the first, toward the [111] axis. The experimental 
results were as follows: in the former structure the plasma 
undergoes Gaussian (symmetric) diffusion, while in the 
latter the diffusion is asymmetric. The authors of Ref. HI 
showed that diffusion data for the tilted heterostructure 
can be fitted by an asymmetric Levy distribution, pre- 
viously obtained^ as a solution to a (one-dimensional) 
generalized Fokker-Planck equation with distinct right- 
and left- diffusion coefficients, and fractional derivatives. 
Furthermore, it was argued that the fractional charac- 
ter of the diffusion is connected to fractal properties of 
the medium, namely a power-law roughness distribution 
of the interfaces that delimit the quantum wells. Fur- 
ther photoluminescence studies suggest that the carrier 
diffusion properties are indeed very sensitive both to in- 
terface roughness, and to the presence of finite- width ter- 
races^.. By varying both optical excitation intensities and 
temperature, a more detailed picture was found for the 
anomalous diffusion phenomena taking place, including 
the likely existence of Auger recombination, which in its 
turn might be induced by fractal interface morphology^. 

The relationship between fractional diffusion and the 
absence of a characteristic length scale (the latter being 
a basic property of fractals) is well-established, and has 
been extensively investigated^i^. 

Here we study a model system along the lines of An- 
derson's picture^i^ d^d^ for electron diffusion in disor- 
dered media. Alternatively to the theoretical approach 



adopted in Refs.QQ, whose starting point was a general- 
ized Fokker-Planck equation for (classical) particles, we 
write a tight-binding Hamiltonian for a (quantum) par- 
ticle, which evolves according to the rules of quantum 
mechanics. The fractal features are incorporated by a 
specific spatial distribution of on-site disorder, to be de- 
scribed below. 

Bearing in mind that the electron-hole pairs undergo- 
ing anomalous diffusion in the experimental systems are 
essentially quantum-mechanical objects, our treatment is 
expected to reflect the basic wave-Hke properties of such 
entities. The link between the quantum properties under- 
lying localization, and those classical features brought to 
fore in classical approaches to diffusion is not well under- 
stood, although it has been discussed in the recent past 
(see, e.g., Ref. [12). 

Our main concern is to provide an independent check 
of whether the connection between fractal properties of 
a disordered medium, and the behavior of particles dif- 
fusing in it, which has been proposed in the classical 
context, is robust enough to survive the translation to 
a quantum-mechanical picture. Therefore, we shall not 
aim at making specific comparisons of our results to the 
experiments described in Refs.[l|I3H. 

In Subsection III A^ we introduce the tight-binding 
hamiltonian and the respective disorder distributions to 
be investigated, along with a brief description of our 
calculational procedures. In Subsection IIIBI results are 
given for the time evolution of the width of wavefunc- 
tion packets against time, as well as for the correspond- 
ing instantaneous profiles. Subsection lll CI deals with the 
introduction of an external bias, in order to investigate 
the effects of anisotropy. Finally, in Sec. Illlj concluding 
remarks are made. 



II. MODEL SYSTEM AND RESULTS 
A. Model 

In order to simplify the calculational framework, we 
restrict our investigation to one-dimensional systems. We 
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consider the one-electron, tight-binding hamiltonian for a 
single-orbital linear chain with nearest neighbor hopping, 



(1) 



where \n) is an orbital with self-energy e„, m — nil, 
and 7 is the hopping energy, assumed to be constant and 
positive. 

We study binary models of disordered alloys, in which 
two different orbitals A (host atoms) and B (substitu- 
tional impurities) coexist, with respective self-energies 
and €3, and overall concentrations 1 — po and po- The 
relevant energy in the present context is e = £5 — 6^1, 
which measures the impurity-host mismatch. 

If the impurity locations are uncorrelated, (a case to 
be denoted as the random binary models or random for 
short) the distances i between any two consecutive im- 
purities have an exponential distribution, i.e.. 



P(£) 



1 -Po 



(l-po)'-exp(-A^) (^=1,2,...)^ 



(2) 

where A — A(po) = — ln(l— Po) is the inverse decay length 
(in lattice parameter units) that sets a scale for the typ- 
ical B ~ B distance, i.e., a cutoff dictating the largest 
allowable distances. Note that the average impurity- 
impurity distance is Pq ^, which approaches A~^ only in 
the low-concentration limit pq where continuum and 
discrete-lattice descriptions become equivalent. For the 
random model, host-impurity {A — B) duality means that 
our investigation can be restricted to < po < 1/2. 

We introduce fractal features by assuming the distri- 
bution of (. to decay with a power law: 



(3) 



where a is a characteristic exponent. 

Distributions similar to Eq. ^ have been exten- 
sively studied in the literature of fractal-based point pro- 
cessesi^iiiii^; in that context, they refer to the probabil- 
ities of occurrence of interevent (time) separations. In 
such cases, most of the interest focuses on the charac- 
teristics of the associated power spectrum, which turns 
out to exhibit 1//" noise properties. Here, by contrast, 
the fractal features of the real-space impurity distribution 
are of interest only inasmuch as they are the background 
against which we simulate the quantum-mechanical evo- 
lution of electronic wavepackets. 

The distribution given by Eq. ^ is normalizable only 
for a > 1; its mean is finite only for a > 2 , and its vari- 
ance is finite only for a > 3. Nevertheless, as discussed 
at length in Subsec. fllB^ it is a physically sensible choice 
to adopt a system-wide normalization which enables one 
to consider 1 < a < 2, for (large but finite) fractal sys- 
tems. In this fractal binary model, or fractal, for short, 
i lacks a typical scale. Here, all distances are allowed, 
on account of the slowly-decaying distribution tail. Thus 
we expect the resulting system to exhibit fractal prop- 
erties, provided that it is large enough. The correlation 



between impurity positions, implied by Eq. ([3]), destroys 
strict A — B duality. 

Comparison between overall concentrations in the two 
cases proceeds by matching the respective average dis- 
tances between impurities, which are p^^ for the former 
model [ as remarked in connection with Eq. ^ ] and 
C(q! — 1)/C(a) for the latter, where C,{a) = Y^^=i " " 
Riemann's zeta function. For example, pq = 0.4 ^ a w 
2.34. Clearly, in the thermodynamic limit any a < 2 
corresponds to vanishing B concentrations in the fractal. 
However, as pointed out above, for large but finite sys- 
tems one can still get nontrivial results for a < 2, once 
suitable normalization considerations are taken into ac- 
count [ see Fig.[3]below, for a full illustration of the po-a 
correspondence, both in the thermodynamic limit and for 
finite systems; see also Eq. |[6]) ]. 



We study the dynamics of wave packet; 
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dom and fractal linear chains. The time evolution of 
the amplitudes V'n(^) — ("-IV'(O) is determined by the 
Schrodinger equation: 



+ 1 



(4) 



where e„ = £11/7 is a dimensionless parameter, and time 
is given in units of ft/7. 

We start with Gaussian wave packets: 



V'„(0) = C cxp 



-{n - n^Y 
2al 



(5) 



where no denotes the initial position of the centroid, and 
(To the initial width, and C is a normalization constant. 

Configurational averages are taken in all the calcula- 
tions presented in this work. 



B. Variances and wave fronts 

We analyze and compare the properties of the two 
kinds of systems just defined. 

We have numerically determined the time evolu- 
tion of the (ensemble-) averaged wave-packets in one- 
dimensional systems of sizes L — 10000, with free ends. 
The number of realizations incorporated in the ensemble 
averages is M = 1000 or, in some cases, larger. 

The calculational method is based on numerical inte- 
gration of Eq. |[4]) with initial condition given by Eq. ^ , 
via a fifth-order Runge-Kutta code. We compared results 
for fixed realizations of disorder, obtained with time steps 
of 0.1 and 0.01, and checked that they are indistinguish- 
able for all practical purposes. Thus, we set the former 
value in all calculations described here. 

The integration is taken up to times not longer than 
enough for the packet to reach the chain's ends, in order 
to avoid refiection effects. The actual hopping-rate of an 
electron in a pure system modeled by Eq. (U is one atom 
per time unit. Considering that diffusion on a disordered 
chain is hampered by impurities, in order to be safe we 
set the upper time limit as L/2, for a packet initially 
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Figure 1: Plots of slopes (time derivatives) d{{Ax)'^)/dt of 
variances of particle probability distribution against time in 
the random system. In (a) e = 1.0 for several po, while in (b), 
po = 0.4 for several e. L = 10 000 and M — 10^ samples each. 



spread over a few sites around the center of an L-atom 
chain. We have used (Tq — 5.0, that is, the initial packets 
are rather locaUzed. 

As usual, we study diflFusional behavior through inves- 
tigation of the time evolution of the ensemble-averaged 
second moment (variance) of the particle probability dis- 
tribution, ((Acc)^). Power-law dependencies, ((Aa;)^ } ~ 
t^^ , with /? = 0, 1/2, and 1 characterize localized, diffu- 
sive, and "ballistic" regimes respectively (the latter to be 
interpreted here in the sense that, even though there is 
no center-of-mass motion in the absence of an external 
bias, the packet width varies linearly in time if each of 
its Fourier components travels unhindered). 

Initially we consider the random model, so as to pro- 
vide a benchmark against which to compare results for 
the fractal case. 

By scanning the values of e = e/7 for the dimension- 
less impurity-host self-energy mismatch, as well as those 
of Po, along suitable intervals, one finds the same overall 
picture, namely apparent ballistic behavior at first, fol- 
lowed by a continuous, smooth crossover towards local- 
ization, i.e. the variances ((Ax)^ ) eventually approach 
saturation. To illustrate this, in Figure [T] we show the 
evolution of the slopes d ((Ax)^ )/c?t against time, in the 
random system. One can see that saturation behavior 
(zero slope) always obtains, albeit at rates which depend 
on £ and po ■ 

We have found no simple scaling picture, from which 
a data collapse plot could be derived. Even for fixed 
e, where for each po one has the characteristic inter- 
impurity distance A~^(po), this quantity does not trans- 
late directly on to a scaling length. In this case, the sat- 
uration value of ((Ax)^ } varies approximately as p^^'^^ 
(for PO < 0.5). 

We now turn to the fractal model, for which represen- 
tative results are shown in Fig. [51 These are for samples 
whose distributions P{t) are normalized with respect to 



Figure 2: Double-logarithmic plots of variance of particle 
probability distribution against time in the fractal system, 
for e = 1.0 and several a (see text for explanation of a < 2 
data). L = 10 000 and M = 5 x 10^ samples each [except 
a = 1.5 for which M — 10^ ]. Inset: asymptotic characteristic 
exponent 2/3 as a function of a, for data on the main diagram 
(see text). 



the actual system size L, i.e., 

^ ' n— 1 

Such a system-wide normaHzation is indeed consistent 
with the inclusion-exclusion principle. This can be seen 
by recalling that, starting from a B atom at x = and 
following along the chain, the statement that, e.g., the 
first B atom occurs at x — £ encompasses all (mutually 
exclusive) 2^~^ possible arrangements of A's and B^s, for 
X > £, and with only A^s for 1 < x < i — 1. With the 
above normalization, even for a < 2 one has, on aver- 
age, a non- vanishing (though in such case rather small) 
impurity concentration. This is illustrated in Fig. [31 

The larger the system, the individual probabilities for 
some given distance £ become comparatively smaller, 
while the occurrence of larger distances becomes possi- 
ble, although with low probability. We find that for larger 
systems, the fractal environment makes it easier for the 
packet to diffuse than the random one, for the same (ef- 
fective) concentration of impurities. Of course, here one 
is only considering diffusion along distances of the order 
of the system's length, L. The crucial difference, relative 
to the random model, is that since there is no intrinsic 
length scale in the fractal system, the single length in- 
troduced by the normalization in Eq. ((Gj) coincides with 
the system size. Thus, although the diffusive behavior 
observed here must always be regarded as an apparent 
regime, it may extend to rather long distances. 

The occurrence of large impurity-impurity distances in 
this system may produce individual realizations of disor- 
der, comprising regions that are effectively pure. The 
result is that the eigenstates in the system may be ex- 
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Figure 3: Effective average concentration po of impurity sites 
in the fractal model, against a. L — > oo refers to Eq. ^ in 
such limit, while points are for L = 10 000, with normalization 
by A'' [Ai' = L is system- wide normalization as in Eq. JG}]; see 
text. 



tended, or semi-extended. Averaging over many realiza- 
tions would produce an ensemble of wave packets com- 
posed of localized and extended states, even though on 
average the system has a finite concentration of impu- 
rities. Least-squares fits of 10^ < i < 5 x 10'^ data of 
Figure [2] to a single power-law form produce an expo- 
nent 2/3 = 1.82(1), 1.30(1), and 1.04(1) respectively for 
a = 1.5, 2.0, and 2.5. Thus one finds an effective anoma- 
lous diffusion regime, which has no counterpart in the 
random disorder model. 

On the other hand, the a = 3.5 results, for exam- 
ple, show that even in the fractal model, one can have 
disorder so strong that locaHzation (within the system's 
limits) is statistically the only possible outcome. For this 
case, the corresponding effective impurity concentration, 
as defined above, is « 0.84, which (for the random model) 
would be equivalent to po — 0.16, by using the po ^ 1—po 
duality valid for the latter type of disorder. One can in- 
fer from Fig. [T] (b) that the behavior of the po = 0.16 
random system is indeed similar to that of the a = 3.5 
fractal case, namely localization setting in for t « 10'^ 
[ we have also checked that the actual saturation value 
of ((Aa:)2), which is « 6 x 10^ for a = 3.5 (fractal), 
falls between the respective ones for po = 0-1 and 0.2 
(random) ]. 

As regards a = 4.0 data, for which one has the equiv- 
alent Po = 0.90 from Fig. [3l comparison is to be made to 
the Po = 0.10 curve of Fig. [T] (b). Again, the agreement 
is very good, as far as general trends are concerned: in- 
creasing a in this range turns out to produce a longer 
localization length. 

The reason why the fractal model behaves similarly to 
the random one for large a, while it certainly does not 
do so for a < 2.5, is as follows. From Eqs. ^ and 
one can work out that the probability distribution for 



distances between consecutive A atoms is 

PA{e) ^ exp[a{a)£], a(a) = lnP(l), (7) 

where P{1) — (j^^{a) is the probability for unit distance 
between B (impurity) atoms. For a < 2.5 (corresponding 
to Po < 0.5), the B atoms are in the minority, so they 
indeed play the role of impurities in an otherwise pure A 
system. As seen above, the fractal properties associated 
to the power-law B — B distance distribution, Eq. 
give rise to the consequent diffusion-like behavior. At 
larger a {po > 0.5), the A atoms are now in the minority. 
Physically, the traveling electrons are only sensitive to 
the existence of two distinct values of on-site energies, 
thus the effective "impurity" label will be assigned in 
practice to the species which occurs less frequently (in 
this case, A atoms, whose distance distribution, Eq. 
is qualitatively the same as in the random model, thus 
bringing about locaHzation). The characteristic length, 
[ln(L(a)]~^, increases with increasing a, thus explaining 
the trend mentioned above. 

In practice, one would expect the region 0.4 < po ^ 
0.6, i.e., 2.35 < a < 2.65, to behave as a crossover re- 
gion. This is because, for the finite systems under study, 
one needs a clear majority of one species over the other 
to be statistically established while still within distances 
shorter than system size. 

Further evidence that effective diffusion-like behavior 
is Hnked to the absence of a typical scale in fractal sys- 
tems can be derived as follows. If, instead of normalizing 
P{i) by system size L as in Eq. we take a fixed N < L 
as the upper Hmit (see in Fig.[3]how this affects the effec- 
tive impurity concentration) , a length scale equal to TV is 
introduced, even though the variation of P{i) against £ is 
still described by a power law. Results for N = 100, and 
assorted values of a are shown in Fig. IH Note that while 
the localization length decreases with increasing a < 2.5, 
the trend is reversed for a > 2.5. This is similar in na- 
ture to the Po ^ 1 — Po duality observed in the random 
model. 

Finally, we examine the actual shapes of particle prob- 
ability densities. For the localized packets in the random 
case, extended exponentials, |-0(a;)P ~ exp(— felxj"^), pro- 
vide reasonably good fits to the region within 20 to 30 
sites from the peak, with 1-5 < (f> < 1.8. At larger dis- 
tances, the probability density decay is somewhat slower 
than that. Fits for 30 < a; < 200 give (j) in the 0.2 - 0.6 
range. 

On the other hand, for the fractal cases with apparent 
diffusion, one can get power-law fits extending to two 
or more decades of distance, as shown in Fig. \E[ Of 
course, the power-law behavior exhibited in the Figure 
is expected to hold only within the system's finite Hm- 
its, as is the case for all features of the apparent diffusion 
regime (recall that our whole study is conducted for times 
not longer than enough for the wavepacket to reach the 
chain's ends, in order to avoid refiection effects). Thus, 
for instance, although data for a = 1.5 (corresponding 
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Figure 4: Asymptotic values of variance of particle proba- 
bility distribution in the fractal system, with P{i) = for 
£ > 100, and the corresponding normalization (see text), for 
£ = 1.0 and several a. L — 10 000 and M — 10^ samples each. 
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Figure 5: Double-logarithmic plots of particle probability 
distribution functions against position, at t = 5000, for as- 
sorted values of a in the (apparent) diffusive regime, e = 1.0, 
L = 10 000, and M — 10^ samples. Lines indicate least- 
squares fits to a power-law form, P{x) ~ . See text for a 
discussion of suitable normalization considerations. 

to power-law decay with S = 0.45) would be strictly non- 
normalizable if such behavior extended to arbitrarily long 
distances, one must keep in mind that (within the present 
context) the wavepacket amplitude will fall to zero before 
reaching the chain's ends. Localized packets for fractal 
disorder, i.e., a > 3.0 (not shown), seem to behave in an 
intermediate way between power-law and extended expo- 
nential decay. 



the Fokker-Planck approach of Refs. UQ. The simplest 
source of such anisotropy is an electric field. 

The effect of applying a uniform electric field to the 
system is incorporated by the inclusion of the term 

Hbias ^ I n)eEn{n\ (8) 

n 

in the Hamiltonian, where e is the electronic charge, and 
E is the field intensity. The corresponding Schrodinger 
equation is 

i^n = (£n+/o")V-'n " 1pn-l " V'n+l , (9) 

where fo = eE/j is a dimensionless bias intensity. 

On a pure system, application of a bias produces a drift 
of the entire wave packet, as its centroid moves through 
the lattice, with diffusion occurring with respect to the 
center-of-mass reference frame. However, due to the un- 
equal effect of the bias on the various Fourier compo- 
nents of the packet, this diffusion may be asymmetric. 
Also, one must be aware of Bloch oscillations, which con- 
fine the wave packet in a region of space, thus producing 
oscillating behavior from application of a static electric 
fieldiS. 

In order to prevent the effect of Bloch oscillations from 
distorting the diffusive behavior which is our main con- 
cern, here we use /o = 1.0 x 10^"^. Then, elementary con- 
siderations show that for a packet starting at the center 
of a chain with L = 10 000, speed reversal will only set 
in at t < 4 X 10'^, giving one a rather broad window of 
observation. 

In a disordered system, the application of a bias gives 
rise to dynamical locaHzation. This is related to Bloch 
oscillations. The difference is that a part of the packet 
diffuses away and performs oscillations, while the other 
part remains localized close to the origin. Such asym- 
metric diffusion behavior may thus be viewed as the co- 
existence of two distinct regimes. 

The results in Fig. [6] show that the diffusing portion of 
the averaged wave packet behaves as if the system were 
pure; its centroid coincides with that of a corresponding 
wave packet in a pure system. One can see that it is only 
for the fractal case with apparent diffusion (a = 2.5) 
that a significant portion of the particle probability 
distribution is pure-system-like. As remarked above, 
for such values of a one has coexistence between indi- 
vidual disorder realizations with effectively delocalized 
eigenstates, and others in which localization occurs. 
The relative height of the corresponding peaks in Fig. [6] 
shows that the latter are much more frequent than the 
former. 



C. Application of a bias 

We search for an effect that, in the present quantum- 
mechanical formulation, would be the equivalent to dis- 
tinct values for left- and right- diffusion coefficients in 



III. DISCUSSION AND CONCLUSIONS 

We have introduced a model system for the incorpo- 
ration of on-site fractal disorder in the one-electron dif- 
fusion problem. As remarked above, the distribution of 
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Figure 6: Particle probability distribution against position, 
for t = 1000, fo = 1.0 X 10-3 for e = 1.0. L = 10 000 and 
M = 10'^ samples each [except a — 2.5 for which M = 5 x 
10^]. Open symbols: fractal systems; full symbols: random 
systems. 



dimensional systems), an observation which is confirmed 
here by the comparison of fractally-disordered and un- 
correlated random electronic systems. 

As regards comparison with experimental results, the 
wavefront profiles depicted in Fig. [6] show a connection 
between anomalous diflFusion of tight-binding electrons 
and fractal properties of the underlying medium. Fur- 
thermore, we have seen that, in the present model, asym- 
metric profiles arise from coexistence of ballistic and lo- 
calized states upon application of an external bias. Such 
coexistence results from the fact that a subset of ensem- 
ble realizations of the (scale-free) disorder are, in fact, 
almost pure. Whether the same explanation holds for 
the experimentally-observed profiles is not certain at the 
moment, though one might conceivably propose ways to 
test it on available samples. 
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